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In this paper we study the spin transport in frustrated antiferromagnetic FCC 
films by Monte Carlo simulation. In the case of Ising spin model, we show that the 
spin resistivity versus temperature exhibits a discontinuity at the phase transition 
temperature: an upward jump or a downward fall, depending on how many parallel 
and antiparallel localized spins interacting with a given itinerant spin. The surface 
effects as well as the difference of two degenerate states on the resistivity are ana- 
lyzed. Comparison with non frustrated antiferromagnets is shown to highlight the 
frustration effect. We also show and discuss the results of the Heisenberg spin model 
on the same lattice. 

PACS numbers: 72.25.-b ; 75.47.-m 

I. INTRODUCTION 

The resistivity in magnetically ordered materials has been studied experimentally and 
theoretically during the last 50 years. Unlike the resistivity in non magnetic systems which 
is due mainly to the scattering of conduction electrons by phonons, the resistivity in mag- 
netic materials depends on the magnetic ordering. It has been theoretically shown that in 
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ferromagnets the resistivity is due to the spin-spin correlation by several workers^"-, using 
various approximations such as mean-field theories and Boltzmann's equation. The main 
difference of these treatments resides on the way the correlation length is taken into account: 
long-range correlation gives rise to a divergence of the magnetic resistivity p at the transition 
temperature Tc, while short-range correlation causes a rounded peak. Experimental data 
show that p has several forms depending on the materials: p shows a peak, a rounded shoul- 
der, or just a change of slope at Tq- In the last case, it is dp/dT which shows a peak^"-. In 
recent experiments performed on different kinds of magnetic pure or doped materials ranging 
from insulators, semiconductors to superconductors, the form of p is very different^"—. These 
theories and experiments suggest that the shape of p depends on the different magnetic in- 
teractions in the system, and on the local magnetic ordering through which the itinerant 
spins evolve. Recently, the study of spin resistivity has attracted again much attention due 
to numerous applications since the discovery of the so-called giant magnetoresistancei^"— . 

Let us recall first some results in the ferromagnetic case. At very low temperatures T 
where spin- wave scattering is dominant, the resistivity has been studied in detail showing 
a behavior— i^^. Near the transition, we have shown in recent papers^^i^^ that the pro- 
nounced peak at the Curie temperature Tc of the magnetic resistivity of ferromagnets can be 
explained by the scattering of itinerant spins by defect clusters formed in the lattice around 
Tc- This picture has an advantage: it can be checked easily by calculating the numbers and 
the sizes of the clusters numerically during the simulation using the Hoshen-Kopelmann's 
algorithn>21. Note that clusters of down spins can be considered as magnetic impurities 
embedded in a up-spin sea which have been theoretically studied by Zarand et al.— where 
they found also a pronounced peak of the resistivity. Of course, the cluster sizes refiect 
the correlation length used in early theories^'^'^. It is not a surprise that our results for 
ferromagnets were in agreement with all these theories. 

The case of antiferromagnets has not been well studied. There were only a few works 
which mentioned briefiy some behaviors. Let us cite the work by Haas where he stated that 
in antiferromagnets there is no peak in p using the spin-spin correlation in the Boltzmann's 
equation'^. Our recent works^^i^ on the simple cubic and the body-centered cubic Ising 
antiferromagnets show that there is indeed no peak in p: p varies with T in a manner 
similar to that of the internal energy versus T. As a consequence, the differential resistivity 
dp/dT shows a peak at the transition temperature Tc just like the specific heat. Interestingly 
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enough, this behavior has been experimentally observed in MnSi and related compounds^i^. 
We will show below some results of these non frustrated antiferromagnets for comparison 
with the frustrated case studied in this paper. 

Our purpose is to show in this paper one of the remarkable cases: the face-centered 
cubic (FCC) antiferromagnet (AF). This system is known to be fully frustrated with a 
strong first-order transition in the Ising case^^i^i. The spin resistivity is known to be very 
sensitive to the nature of the ordering of the media through which the itinerant spins move: 
local disordering (for instance, disordering near film surfaces, around magnetic impurities), 
magnetic instability, ... The FCC AF is thus a very good candidate where exotic behaviors 
are expected for the spin resistivity. This will indeed be seen in this work. By using Monte 
Carlo (MC) simulation, we show that frustration and interaction range between itinerant 
spins and lattice spins play a crucial role on the spectacular discontinuity of the magnetic 
resistivity at Tc- In the case of Heisenberg spin, the transition is also of first order, though 
weaker, in the bulk—. We study also this model in this work to outline the effect of spin 
continuous degrees of freedom on the resistivity. 

The paper is organized as follows. Section II is devoted to the description of the model and 
the calculation method. In section III, we show MC results on the temperature dependence 
of the magnetic resistivity in the Ising spin model. Section IV is devoted to results of 
the Heisenberg case. Discussion and explanation are given with regard to the transport 
mechanism. Concluding remarks are given in section IV. 

II. MODEL AND SIMULATION PROCEDURE 

A. Model 

We consider a thin film of FCC lattice structure where each lattice site is occupied by 
an Ising spin whose values are ±1. The Heisenberg spin model is considered in section 
IIVI Interaction between the lattice spins is limited to nearest-neighbor (NN) pairs with the 
following Hamiltonian : 

V-i = JijSi.Sj (1) 
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where Si is an Ising spin, Jjj the exchange integral between the NN spin pair Si and Sj. 
Hereafter we take Jij = Jg for surface spins and Jjj- = J for other NN spin pairs. We 
consider here the antiferromagnetic interaction J^^J < for the rest of this paper. The 
system size is x Ny x where A^^ is the number of FCC cells in the x direction etc. 
Periodic boundary conditions (PBC) are used in the x and y directions while the surfaces 
perpendicular to the z axis are free. The film thickness is A^^. 

The FCC AF is a fully frustrated system which is composed of tetrahedra each of which 
has four equilateral triangles. We know that it is impossible to fully satisfy simultaneously 
the three antiferromagnetic bond interactions on each triangle. As a consequence, the bulk 
lattice has an infinite ground-state degeneracy- . In the case of a thin film, the surface spin 
configuration depends on Jg as shown in Fig. O^. 




FIG. 1: Ground state spin configuration of the FCC cell at the film surface (basal xy plane). The 
horizontal (vertical) axis is the x (z) axis. Upper: ground state when \ Js\ < 0.5|J|, Middle and 
Lower: first and second ground states when | Js| > 0.5| J|. 



For < 0.5 1 J|, the ground state is composed of ferromagnetic xy planes antiferromag- 
netically stacked in the z direction a shown in the upper figure of Fig. [T] For \ Js\ > 0.5| J|, 
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the ground state is two-fold degeneracy as shown in the middle and lower figures of Fig. [H 
The difference of these two configurations is that the middle figure is an alternate stacking 
of up- and down-spin planes in the y direction while the lower figure is an alternate stacking 
of up- and down-spin planes in the x direction. These degenerate states are not equivalent 
in the spin transport in the x direction as seen below: in the first degenerate state, the 
itinerant spins move in the x direction between an up-spin plane and a down-spin plane, 
while in the second degenerate state the itinerant spins meet successively an up-spin plane 
and a down-spin plane perpendicular to their trajectories. We will present our results for 
these two cases separately. 



B. Multi-step Averaging 

The procedure of our simulation can be split into two steps. The first step consists 
in equilibrating the lattice at a given temperature T without itinerant electrons. When 
equilibrium is reached, we study thermodynamic properties of the film so as to determine its 
Neel temperature by examining quantities like internal energy, specific heat, susceptibility 
and magnetization as functions of T— i^I. 

In the second step, we randomly add A^o polarized itinerant spins into the thermalized 
FCC lattice. In the structure, each itinerant electron interacts with lattice spins in a sphere 
of radius Di centered on its position, and with other itinerant electrons in a sphere of radius 
D2- We define these interactions as follows 

Hr = -Y.kA-S, (2) 

where cXi is the Ising spin of itinerant electron and Jj.-,- denotes the interaction that depends 
on the distance between an electron i and the spin Sj at the lattice site j. We use the 
following interaction expression : 

lij = V-"^'^ with rij = \fi-fj\ (3) 

where /q and a are constants which will be chosen in section III CI In the same way, inter- 
action between itinerant electrons is defined by : 

'Hm = -^Kijdi.dj (4) 

K,^, = K^e-^'^^ (5) 
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with (Tj the spin of itinerant electron and Kij the interaction that depends on the distance 
between electrons i and j. The choice of the constants Kq and /3 is discussed in III CI 
Dynamics of itinerant electrons is ensured by an electric field applied along the x axis. 
Electrons enter the system at the first end, travel in the x direction, leave the system at the 
second end. The PBC on the xy planes ensure that the electrons who leave the system at 
the second end are to be reinserted at the first end. For the z direction, we use the mirror 
reflection at the two surfaces. These boundary conditions are used in order to conserve the 
average density of itinerant electrons. One has 



where e is the charge of electron, e the applied electrical field and ^ the displacement vector 
of an electron. 

Since the interaction between itinerant electron spins is attractive, we need to add a chemical 
potential in order to avoid a possible collapse of electrons into some points in the crystal 
and to ensure a homogeneous distribution of electrons during the simulation. The chemical 
potential term is given by 



where n(r) is the concentration of itinerant spins in the sphere of D2 radius, centered at r. 
D is a. constant parameter appropriately chosen. 

The procedure of spin dynamics can be described as follows. After injecting A^o itinerant 
electrons in the equilibrated antiferromagnetic FCC lattice, we calculate the energy Eoid of 
an itinerant electron taking into account all interactions described above. Then we perform 
a trial move of length I taken in an arbitrary direction with random modulus in the interval 
[Ri^R2\ where -Ri = and R2 = a/\/2 (nearest-neighbor distance), a being the lattice 
constant. Note that the move is rejected if the electron falls in a sphere of radius tq centered 
at a lattice spin or at another itinerant electron. That excluded space emulates the Pauli 
exclusion. We calculate the new energy Emw and use the Metropolis algorithm to accept 
or reject the electron displacement. We choose another itinerant electron and begin again 
this procedure. When all itinerant electrons are considered, we say that we have made a 
MC sweeping, or one MC step. We have to repeat a large number of MC steps to reach 



= -eel 



(6) 



Uc = DVrn{f) 



(7) 
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a stationary transport regime. When the stationary regime is reached, we perform the 
averaging to determine physical properties such as magnetic resistivity, electron velocity, 
energy etc. as functions of temperature. 

We emphasize here that in order to have sufficient statistical averages on microscopic 
states of both the lattice spins and the itinerant spins, we use the following procedure: after 
averaging the resistivity over Ni MC steps we thermalize again the lattice with N2 steps, 
then equilibrate the itinerant spins with A^3 steps before taking back the averaging of the 
resistivity during Ni steps. We repeat this cycle {Ni + A^2 + ^3) for A^4 times (typically 
iV4=100). The total MC steps is therefore equal to N4 x (A^i + A^2 + N^). As will be 
seen below, this procedure reduces strongly thermal fluctuations observed in our previous 
work^-. The transport averaging is made for N2 = 100 configurations of lattice spins. At 
each configuration, we adjust A^^i so that the spin resistivity is calculated during 1000 lattice 
sweepings per itinerant spin (each electron passes through the system 1000 times). N2 
depends on T: it can be several thousands near Tc- In all, at each T the initial equilibration 
time for lattice spins lies around 10^-10® steps per spin and statistical averages are made 
with about 4 x 10^ steps per spin (A^4 x Ni). 

We define resistivity p as : 

P = - (8) 

He 

where rig is the number of itinerant electrons crossing a unit slice perpendicular to the x 
direction per unit of time. 

In this paper we use the lattice size = Ny = 20 and A^^ = 8. 

For studying the spin transport, we consider A"o = (A^^ x Ny x Nz)/2 itinerant spins 
(one electron per two FCC unit cells). Except otherwise stated, we choose interactions 
Iq = Ko = 0.5, Di e [0.6a; 2a], D2 = a, D = 0.35, e = 1, A'o = 1600, and tq = 0.05a. A 
discussion on the effect of a variation of each of these parameters will be given in III CI 

Note however that, due to the form of the interaction given by Eq. the itinerant 
spins have a tendency to form compact clusters to gain energy. This tendency is neutralized 
more or less by the concentration gradient term, or chemical potential, given by Eq. ([7]). 
The value of D has to be chosen so as to avoid a collapse of itinerant spins. We show in Fig. 
|2]the phase diagram in the space {Kq, D). The limit depends of course on the values of Di 
and D2. 
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C. Choice of different parameters 

We will show below results obtained for typical values of parameters. The choice of 
the parameters has been made after numerous test runs. Let us describe the principal 
requirements for the choice: 

i) we choose the interaction between lattice spins as unity, i. e. |J| = 1. 

ii) we choose interaction between an itinerant and its surrounding lattice spins so as its 
energy Ei in the low T region is the same order of magnitude with that between lattice 
spins. To simplify, we take a = 1. 

iii) interaction between itinerant spins is chosen so that this contribution to the itinerant 
spin energy is smaller than Ei in order to highlight to effect of lattice ordering on the spin 
current. To simplify, we take /3 = 1. 

iv) the choice of D is made in such a way to avoid the formation of clusters of itinerant 
spins (collapse) due to their attractive interaction [Eq. (^] as shown above. 

v) the electric field is chosen not so strong in order to avoid its dominant effect that would 
mask the effects of thermal fluctuations and of the magnetic ordering. 

vi) the density of the itinerant spins is chosen in a way that the contribution of interactions 
between themselves is neither so weak nor so strong with respect to Ei. 

Within these requirements, a variation of each parameter does not change qualitatively 
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the results shown below. As will be seen, only the variation of Di does change drastically 
the results. That is the reason why we will study in detail the effect of this parameter. For 
larger densities of itinerant spins, the resistivity is larger as expected because of additional 
scattering process between itinerant spins. 

In view of the above requirements, we take for the simulations: J = — 1 (AF interaction), 
Iq = Ko = 0.5|J|, Di e [0.6a; 2a], D2 = a, D = 0.35, e = I, Nq = 1600, and tq = 0.05a, 
a being the FCC lattice constant. Within these choices, the results in the following will be 
presented in the following units: the spin energy is in the unit of \J\, the temperature is in 
the unit of \J\/kB, the distance is in the unit of a. 

Finally, we keep a constant when varying Di. This is because varying Di means we 
include or not include some far neighbors. If the interaction of these far neighbors follows 
the same Eq. (3) as the one between shorter neighbors then it is known theoretically 
that no new interesting physical can occur except the modification of non universal values 
such as the critical temperature. This case corresponds to ferromagnetic interaction where 
further neighbors do not cause interesting effect. On the other hand, when further neighbor 
interaction is in competition with the interaction of nearer neighbors, the system can be 
frustrated, then physical properties can radically vary. This is the case studied here. Now 
varying a to keep the the interaction constant will not change observed physical behaviors. 
What will change is the relative value of the energy and therefore the value of the transition 
temperature but not the qualitative behavior of the system. 

III. RESULTS FOR THE ISING CASE 

We show in Fig. [3] the staggered magnetization of the lattice as a function of T. As seen 
here the transition is of first order with a discontinuity at Tc — 1.79. Note that the Ising 
AF FCC thin film shows a first-order transition down to a thickness of about four atomic 
layer s2i. 

A. Resistivity in the first degenerate state 

We consider the first degenerate configuration shown in the middle figure of Fig. [1] with 
Jg = J = —1. To understand the behavior of the spin resistivity which will be shown below. 
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FIG. 3: Staggered magnetization of antiferromagnetic FCC thin film of thickness A^^ = 8 versus 
T. The transition temperature Tc — 1.79. 

let us first show how various physical quantities at a given temperature depends on Di. 

In the ferromagnetic state, increasing (decreasing) Di results in an increase (decrease) 
of the number of "parallel" lattice spins which interact with the itinerant spin. This means 
that increasing (decreasing) Di results in a decrease (increase) of the energy of the itinerant 
spin. In antiferro magnets, the situation is different: changing Di will result in an oscillatory 
change of the difference of the numbers of parallel and antiparallel spins in the sphere of 
radius Di, namely AN^^ = ~ N^. This is because of antiferromagnetic ordering. 

We show in Fig. Hlthe resistivity, the spin velocity in the x direction, AA^^^ and the energy 
of an itinerant spin at T = 1.65 below the transition and at T = 2 in the paramagnetic 
lattice, for different values of Di. The following remarks are in order: 

• As said above, at low T, AA^-^^ oscillates with varying Di. When AN-^^ is maximum, 
i. e. the number of up spins is large, the energy of the itinerant spin is low. As a 
consequence the itinerant spin will not move easily under the electric field, its velocity 
is therefore slowed down, making the resistivity to increase. 

• At high T, the lattice spins are disordered, there are no more shells alternately of up 
spins and down spins around an itinerant spin. So the oscillatory behavior of AA^^j^ is 
reduced as seen in Fig. |l]at T = 2. 

In Fig. Owe show the spin resistivity p versus T for two typical values of Di. In all 
cases resistivity p is small for low T then increases with increasing T. At T^, it undergoes 
a discontinuity upward jump. After transition, the resistivity decreases slowly to the same 
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FIG. 4: Different pliysical quantities versus Di in unit of the lattice constant a in the case of the 
first degenerate configuration. From top to bottom: Resistivity, velocity on the x axis, difference 
of up- and down-spin numbers, energy of an itinerant spin. In each plot circles corresponds to 
T = 1.65, and diamonds to T = 2. A^^ = 8, A^o = 1600, Js = J = -1.0 and D = 0.35. 



value for all Di in paramagnetic phase. We explain the behavior of p at different temperature 
regions: 

• When T — 7- 0, the resistivity slightly increases because itinerant spins search to mini- 
mize energy by occupying low-energy positions in the periodic lattice. Since thermal 
energy and electric field are not strong enough to make them move, the itinerant spins 
are somewhat frozen in some almost periodic positions, namely a pseudo crystalliza- 
tion. We have studied the spatial distribution of itinerant spins. The results show 
indeed a radial distribution with peaks up to rather long-range positions at low T, 
namely up to 4th nearest neighbors. Note that the increase of p when T — )■ has been 
observed in many experiments among which we can mention: Fig. 11 of the paper by 
Chandra et al. on CdMnTe,— Fig. 2 of the paper by Du et al. for MnFeGe,— Fig. 
6a of the paper by McGuire et al. on AF superconductors LaFeAsO,^^ Fig. 2 of the 
paper by Lu et al. on AF LaCaMnO,-^^ and Fig. 7 of the paper by Santos et al. on AF 
LaSrMnO.^° Note however that most of these experiments concern doped semiconduc- 
tors. In semiconductors, the carrier concentration increases with increasing T. Our 
model has a number of itinerant spins which is independent of T in each simulation. 
So we cannot compare quantitatively our results with experiments on semiconductors, 
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in particular these latter are often magnetically disordered systems. But in each simu- 
lation, we can take another concentration (see our previous paper—): the results show 
that the resistivity is somewhat modified but keeps the same feature, except the fact 
that the stronger the concentration is the smaller the peak at Tc becomes. Therefore, 
we believe that some generic effects independent of carrier concentration will remain. 

Note also that our results of the resistivity at low T depends on our model. This 
behavior (p increases with decreasing T) is observed only when we introduce a rather 
strong interaction between itinerant spins (variable Kq). Reducing Kq will suppress 
this tendency. 

On the hypothesis of frozen electrons, there is a reference on the charge-ordering at 
low T in Pro.sCao.sMnO^ due to some strain interaction. A magnetic field can make 
this ordering melted giving rise to a depressed resistivity. Though our model does 
not correspond to this material, the fundamental concept is similar. For the system 
Pro.sCao.sMnOs, which shows commensurate charge order, the "melting" fields at low 
temperatures are high, on the order of 25 Teslai^. 

We mention here that low-T behaviors can be also studied by alternative Kubo and 
Landauer methods as shown by Ref. 22. 

When T increases, thermal energy unfreezes itinerant electrons, the system is progres- 
sively unfrozen and the resistivity slightly decreases and then increases up to Tc- 

• At Tc, p exhibits a discontinuity due to the discontinuity of the lattice magnetization 
to which the itinerant spins are coupled. For the first degenerate configuration (Fig. 
dl middle) p makes an upward jump. 

• After Tc, the lattice is paramagnetic: there is no significant effect of Di as discussed 
earlier. 

We show now the effect of the magnetic field in Fig. O We observe here that the peak 
hight increases with increasing B, contrary to the case of ferromagnets where the peak 
diminishes with increasing B^^. The difference can be explained by the fact that in anti- 
ferromagnets the magnetic field causes a transition at some T by returning antiparallel spins 
and thus enhances critical fiuctuations while in ferromagnets the magnetic field suppresses 
fiuctuations and forbids a phase transition. Since the peak height is proportional to critical 
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FIG. 5: Resistivity of thin film of size Nx = Ny = 20 and A''^ = 8 for A'^o = 1600 itinerant spins 
versus T for Di = a (black circles) and Di = 1.25a (white circles), a being the lattice constant. 
Case of the first degenerate state. Jg = J = —1.0, Iq = Kq = 0.5, D = 0.35. 



fluctuations, it is not surprising that the peak increases with increasing B in antiferromag- 
nets. Note that Tc diminishes with increasing B as expected in antiferromagnets. 
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FIG. 6: Resistivity versus T for two values of magnetic field B with Di = a, = 8, A'^o = 1600 
and Js = J = — 1.0. Black circles correspond to i? = 0.75 and white circles to i? = 0.25. 
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B. Resistivity in the second degenerate state 

Let us consider the second degenerate configuration where the ferromagnetic up- and 
down-spin planes are perpendicular to the spin current in the x direction (lower figure in 
Fig. [1]). We show first in Fig. [7] the resistivity, the electron velocity, AA''^^ = — Ni and 
the energy of an itinerant spin at two temperatures, below and above Tc, as functions of 
Di. One observes here a cross-over of low-T and high-T resistivities at different positions 
of Di. So, depending on Di, low-T resistivity can be smaller or larger than that of high-T. 
At Tc, p can jump upward or downward depending on the value of Di. This is what we 
see in Fig. [HI Note that, as in the case of the first degenerate configuration, a minimum 
energy of itinerant spin corresponds to a minimum of the velocity and a maximum of the 
resistivity. These quantities are closely related to each other as expected from the physical 
picture described above. 
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FIG. 7: Resistivity, spin velocity, AA^^ and energy of itinerant spin versus Di (in unit of the lattice 
constant a) at temperatures T = 1.65 (black circles) and T = 2.0 (white diamonds) in the case of 
second degenerate configuration. Nz = 8, Nq = 1600, Jg = J = —1.0. 



We can approximately identify the intervals of Di where p jumps (falls) by looking at 
the top panel of Fig. [T) at a given Di, p jumps (falls) when p at T = 1.65 is lower (higher) 
than that at T = 2. 
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FIG. 8: Resistivity versus temperature in the case of second degenerate state for Di = a (black 
circles) and Di = 1.25a (white circles) with = S, Nq = 1600, = J = -1.0, Iq = Kq = 0.5, 
D = 0.35. 
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FIG. 9: Resistivity versus T for first (black points) and second degenerate (white points) configu- 
rations for Di = a with A^^ = 8, iVo = 1600, Js = J = -1.0, lo = Kq = 0.5, D = 0.35. 



C. Surface Effects 



In order to enhance the surface effect, in addition to a small value of Jg we allow the 
exchange interaction between a surface spin and its neighbors in the beneath layer to be Jp 
which will be taken to be small in magnitude. We show in Fig. [TD]the surface magnetization 
and the magnetizations of the interior layers as functions of T for Jg = Jp = —0.5 and 
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J = — 1. As seen here, the surface transition takes place at a lower temperature Ti ~ 1.2 
while interior layers become disordered at T2 ^ 1.8. As a consequence, one expects that 
the surface fluctuations at Ti will induce an anomaly in p in addition to that at T2. This is 
shown in Fig. [TTJ Note that the increase of p at low T is an effect of a pseudo crystallization 
of itinerant spins at low T as discussed above. 
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FIG. 10: Layer magnetizations versus T for Jg = Jp = —0.5 and J = —1. Other parameters: 
Di = a, = 8, iVo = 1600, Iq = = 0.5, D = 0.35. The surface transition is at Ti ~ 1.2. 
The vertical dotted hne is a guide to the eye indicating the discontinuous fah of interior layer 
magnetization. 




D. Traveling paths 

Let us show now how the itinerant spins choose their paths to travel across the lattice. 
We show in Fig. [12] the energy landscape at T = 1 for both degenerate configurations. This 
gives some information concerning the spatial energy distribution in the system. 

As we said above, the spin motion depends solely on the spin energy due to its interaction 
with surrounding spins. The lower energy it has the longer it stays in that position. By 
examining the different traveling paths we come to this observation: for a given Di, the 
itinerant spin will choose its path where its energy is low. This is understandable from a 
viewpoint of statistical physics. Paradoxically, by choosing low-energy paths, its motion is 
slowed down because as said earlier itinerant spins feel energetically at ease so it does not 
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FIG. 11: Resistivity versus temperature T in the case shown in Fig. [TOl There are two anomahes 
occurring respectively at the surface transition temperature and at the bulk one. 




FIG. 12: Energy landscape at T = 1 in a cubic box of 2a x 2a x 2a dimension where a is the lattice 
parameter, for the first and second degenerate configurations (left and right, respectively). The 
energy scale is indicated on the figure. 



want to move. So, depending on the value of Di, itinerant spins will move near up-spin 
planes or near down-spin planes in order to have a low energy (see Fig. [T^ . 

Let us show now in Fig. [T3] how the spins travel across the system. As seen, for equal 
travel time an itinerant spin moves faster and farther in the first degenerate configuration 
than in the second one for Di = a. This can be understood because the itinerant (up) spin 
is stopped for a more or less long time in front of a wall of down spins perpendicular to 
its X trajectory in the second degenerate configuration. However, when Di is very large, 
for instance 1.4a, there is no more difference between the two configurations because the 
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FIG. 13: (Color on line) 3D antiferromagnetic FCC lattice at T = 1, where lattice down spins are 
presented in red, up spins in yellow and itinerant spins (which are up) in white. The plane shown 
is the xy plane with x direction (spin flow direction) being horizontal. Left: snapshot in the case 
Di = a. Right: snapshot in the case Di = 1.4a. 

distance is long enough for an itinerant spin to see other spins across down-spin walls. 

E. Discussion 

Let us recall that in ferromagnets, p shows a peak at Tc- This peak was interpreted as 
a consequence of spin-spin correlation. The form of the peak depends on the correlation 
range^'^'^. Other interpretations which were based on scattering by defect clusters^^ or by 
impuritiea^^ are all in agreement. In antiferromagnets, the situation is quite different. Unlike 
in ferromagnets where itinerant spins are slowed down only when they encounter antiparallel 
spins of defect clusters, in antiferromagnets itinerant spins see both parallel and antiparallel 
spins in any of its position, so their motion depends drastically on their immediate local 
spin configuration whose energy is determined by the interaction range Di. In addition, 
the behavior of the spin resistivity in antiferromagnets depends on several other ingredients 
among which one can mention the crystal structure, the nature of the magnetic ordering, 
and the spin model. We have simulated some non frustrated antiferromagnets such as Ising 
antiferromagnetic simple cubic (SC) and body-centered cubic (BCC) lattices. The spin 
resistivity shows no peak in these cases^^. 
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FIG. 14: Travel path of an itinerant spin at T = 1 in the first degenerate state (upper) and in 
the second degenerate state (lower) during an equal travel time with Di = a. As seen, spins move 
more easily in the first degenerate configuration than in the second one. Other parameters are 
iVo = 1600, Iq = Ko = 0.5, D = 0.35, Js = J = -1.0 
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IV. RESULTS FOR THE HEISENBERG CASE 

In this section, we presently briefly the results on the same lattice with the Heisenberg 
spin model. Itinerant spins are the same as used above, namely polarized Ising spins. This 
assumption allows to outline only the effect of the continuous nature of the Heisenberg lattice 
spin on the resistivity. The full Hamiltonian with different kinds of interaction is assumed 
as above except the exchange interaction between lattice spins. This is given by 

'H = -J2J^,S,-S,-AJ2S!S^ (9) 

where Sj is the Heisenberg spin at the site i and A an Ising-like anisotropy which is assumed 
to be negative to favor an antiparallel spin ordering on the z axis. When A is zero, one has 
the isotropic Heisenberg model. In order to have at phase transition at a nonzero T, we 
should take a nonzero value for A because it is known, by the theorem of Mermin- Wagner— , 
that for vector spin models there is no long-ranged ordering at finite temperatures in two 
dimensions. The small thickness considered here is, in a phase-transition point of view, 
equivalent to a two dimensional system. Except A, note that we use the same assumptions 
as in Eq. ([1]). 

The transition temperature with A = —1 is T^^ ~ 0.79 for the lattice size = Ny = 20, 
Nz = 8. We use here the same analysis as for the Ising case above: we first look at the 
effect of Di on the resistivity at two temperatures, one lower and one higher than Tc- This 
is shown in Fig. [15] where the upper (lower) figure is for the first (second) degenerate spin 
configuration. Again here, one observes that the two degenerate states do not yield the same 
transport properties as in the Ising case. The same remarks on physical mechanism are thus 
applied fsee IIII Al and IIIIAI) . 

Let us show now in Fig. [12] the resistivity as a function of T for two typical values of 
Di. As seen, depending on the value of Di, p undergoes a sharp increase or decrease at 
Tc- At some values such as that corresponding to the upper curve of the upper figure, 
the resistivity can go across a large region of fluctuations without a sharp jump. So in 
experiments, care should be taken to interpret similar behavior if any. Note that the second 
degenerate configuration yields always a larger resistivity than in the first one, as observed 
in the Ising case in the previous section. 

The effect of A on the resistivity is not very important in the reasonable range [0.1, 1.5]: 
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FIG. 15: Heisenberg case. Resistivity, spin velocity, and energy of itinerant spin versus Di (in unit 
of the lattice constant a)at temperatures T = 0.75 (black circles) and T = 0.85 (white diamonds) 
for first (upper) and second (lower) degenerate configurations. A = —1, A^^ = 8, Nq = 1600, 
J, = J = -1.0. 



except the fact that Tc varies with A, for instance Tc — 0.65 for A = 0.5 and Tc — 0.55 for 
A = 0.1, the discontinuity of p at Tc diminishes only slightly with decreasing A. 

V. CONCLUSION 

The model used in this paper is rather general. It has been applied to ferromagnets with 
success^^. We believe that it can be applied to different materials by choosing appropriate 
interactions. For example, in metals we have to reduce to almost zero the interaction between 
itinerant spins and lattice spins to create the situation of almost-free conduction electrons. 
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FIG. 16: Heisenberg case. Resistivity of thin film of size = Ny = 20 and A^^ = 8 for A'^o 
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itinerant spins versus T for Di = a (black circles) and Di = 1.25a (white circles) in unit of the 
lattice constant a for first (upper) and second (lower) degenerate states. A = —1, Jg = J = —1.0, 
Iq = Ko = 0.5, D = 0.35. 



The lattice disordering transition in magnetic metals then should not strongly affect the 
spin resistivity. In semiconductors, that interaction should be strong enough to reproduce 
a peak of p as experimentally observed. 

We have shown in this paper that the spin resistivity p of the fully frustrated 
FCC antiferromagnet is quite different from that of ferromagnets^^ and non frustrated 
antiferromagnets^. p does not show a peak at the magnetic phase transition temperature. 
It shows instead a discontinuous jump at the transition temperature Tc- The jump depends 
on the numbers of parallel and antiparallel localized spins which interact with an itinerant 
spin. After transition, the resistivity tends to a saturation value independent of Di. The 
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abrupt behavior of p at Tc in the AF FCC Ising lattice is an effect of the frustration which 
causes a first-order transition of the lattice magnetic ordering leading to a discontinuity of 
p at Tc- 

We are not aware of experiments performed on spin transport in materials with first- 
order magnetic transition. Our result is thus a prediction which would be useful for future 
experiments. Note however that for electrical transport, the electrical resistivity shows a 
discontinuity at a metal- insulator "first-order" transition in PrNiO^^ and NdNiO^^. Our 
magnetic resistivity found in this paper has also a discontinuity behavior at a magnetic 
"first-order" transition. This similarity shows that the resistivity is closely related to the 
nature of the phase transition, whatever its origin (magnetic, insulator- metal, ...) may be. 
The mapping between the two cases however is not the scope of this paper. 

We have also shown that the surface disordering causes a peak of the resistivity at the 
surface transition temperature. In the Heisenberg model, the spin continuous degrees of 
freedom weaken the first-order transition, yielding in general a reduction of the critical 
temperature and a less abrupt change of the resistivity at the transition. 

As a last remark, let us emphasize that the behavior of the spin resistivity at Tq is quite 
different from one antiferromagnet to another. It depends on many factors such as the 
lattice structure, the interaction range, the spin model and the instability (in particular due 
to frustration) of the spin ordering. We have studied here the effects of some of them, but 
a throughout understanding needs much more investigations and analysis. 
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